nature  v>= — 

COMMUNICATIONS 


ARTICLE 


Received  16  Sep  2015  |  Accepted  25  May  2016  |  Published  22  Jun  2016 


DOI:  10.1038/ncomms12045 


OPEN 


Experimental  demonstration  of  the  microscopic 
origin  of  circular  dichroism  in  two-dimensional 
metamaterials 


A.B.  Khanikaev1,  N.  Arju2,  Z.  Fan2,  D.  Purtseladze2,  F.  Lu3,  J.  Lee3,  P.  Sarriugarte4,  M.  Schnell4,  R.  Hil lenbrand4, 
M.A.  Belkin3  &  G.  Shvets2 


Optical  activity  and  circular  dichroism  are  fascinating  physical  phenomena  originating  from 
the  interaction  of  light  with  chiral  molecules  or  other  nano  objects  lacking  mirror  symmetries 
in  three-dimensional  (3D)  space.  While  chiral  optical  properties  are  weak  in  most  of  naturally 
occurring  materials,  they  can  be  engineered  and  significantly  enhanced  in  synthetic  optical 
media  known  as  chiral  metamaterials,  where  the  spatial  symmetry  of  their  building  blocks  is 
broken  on  a  nanoscale.  Although  originally  discovered  in  3D  structures,  circular  dichroism  can 
also  emerge  in  a  two-dimensional  (2D)  metasurface.  The  origin  of  the  resulting  circular 
dichroism  is  rather  subtle,  and  is  related  to  non-radiative  (Ohmic)  dissipation  of  the 
constituent  metamolecules.  Because  such  dissipation  occurs  on  a  nanoscale,  this  effect  has 
never  been  experimentally  probed  and  visualized.  Using  a  suite  of  recently  developed 
nanoscale-measurement  tools,  we  establish  that  the  circular  dichroism  in  a  nanostructured 
metasurface  occurs  due  to  handedness-dependent  Ohmic  heating. 
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Just  as  left  or  right  hands  that  lend  the  name  to  the  fascinating 
geometric  concept  of  chirality,  certain  types  of  molecules 
cannot  be  brought  into  congruence  with  their  mirror  image 
enantiomer.  From  quartz  crystals  to  organic  sugars  and  nucleic 
acids,  such  ‘chiral’  molecules  are  ubiquitous  in  the  physical  world. 
For  reasons  not  yet  entirely  understood,  biological  molecules  tend 
to  occur,  predominantly,  in  either  left-  or  right-handed  chiral 
forms1.  Similarly,  due  to  their  vector  nature,  light  waves  also 
naturally  possesses  two  forms  of  handedness  and  can  exist  as  left 
and  right  circularly  polarized  (LCP  and  RCP,  respectively)  waves. 
So  it  is  not  surprising  that  one  of  the  earliest  discovered 
manifestations  of  chirality  belongs  to  optics:  chiral  materials  are 
responsible  for  the  phenomenon  of  optical  activity  that  results  in 
the  rotation  of  the  polarization  state  of  light.  Polarization  rotation 
occurs  because  light  of  opposite  handedness  perceives  different 
refractive  indices  and  absorption  coefficients  when  propagating 
through  any  homochiral  material.  These  effects  manifest 
themselves  as  unequal  absorption  rates  (circular  dichroism) 
and  unequal  phase  advances  (circular  birefringence)  for  LCP  and 
RCP  light. 

Because  chiral  response  in  naturally  occurring  materials  is 
quite  weak,  with  circular  birefringence  of  the  order  of  fine 
structure  constant  a  =  1/137  (refs  2,3),  alternative  approaches 
based  on  metamaterials  and  metasurfaces  have  been  recently 
explored4.  The  most  straightforward  way  of  producing  strong 
circular  dichroism  and  birefringence  is  to  mimic  natural  helical 
molecules,  such  as  glucose,  by  synthesizing  fully  three- 
dimensional  (3D)  chiral  metamolecules  that  look  like  helices5-8 
Such  artificial  structures  exhibit  strong  transmission  difference 
between  LCP  and  RCP  polarization  states  of  light,  as  well  as 
handedness-dependent  non-linear  response9.  However,  an 
alternative  approach  to  achieving  circular  dichroism  that  does 
not  rely  on  3D  nature  of  the  metamolecules  has  recently  emerged. 
Fully  two-dimensional  (2D)  metamaterials,  also  known  as  meta¬ 
surfaces10-18,  comprised  of  planar- chiral  plasmonic  meta¬ 
molecules  that  are  just  nanometres  thick,  have  been  shown  to 
exhibit  chiral  dichroism  in  transmission  (CDT).  Theoretical 
calculations 14,1 9,20  indicate  that  this  surprising  effect  relies  on 
finite  non -radiative  (Ohmic)  losses  of  the  metasurface.  In  the 
absence  of  such  losses  on  the  nanoscale,  the  CDT  defined  as  the 
difference  A  TCD  =  TR—  TL  between  the  transmission  coefficients 
Tr  and  Tl  of  the  RCP  and  LCP  waves  must  identically  vanish. 
This  surprising  theoretical  prediction  has  never  been 
experimentally  verified  because  of  the  challenge  of  measuring 
non-radiative  loss  on  the  nanoscale. 

In  this  study  we  use  a  combination  of  nanoscale- characteriza¬ 
tion  techniques  to  demonstrate  that  the  RCP  and  LCP  states  of 
the  incident  light  produce  drastically  different  distributions  of 
optical  energy  and  Ohmic  heat  dissipation  in  the  2D  chiral 
nanoantennas,  thereby  producing  a  strong  chiral  dichroism  in 
absorption  (CD A)  AACD  =  AR  —  AL.  A  planar-chiral  metasurface, 
along  with  its  chiral  enantiomer,  was  designed  to  maximize  the 
CD  A  in  mid-infrared  wavelength  range21.  The  CD  A  gives  rise  to 
the  CDT  observed  experimentally  in  the  far-held  measurements. 
We  then  use  scattering- type  near-field-scanning  optical 
microscopy  (s-SNOM)22,23  to  map  the  optical  energy 
distribution  on  the  nanoantennas  and  their  enantiomers  in 
response  to  the  RCP  and  LCP  light.  Photo-expansion 
microscopy24-28,  also  known  as  AFM-IR,  was  then  utilized  to 
experimentally  demonstrate  drastically  different  Ohmic  heating 
of  the  nanoantennas  under  RCP  and  LCP  light  illumination. 


Results 

Design  and  theoretical  analysis  of  the  metasurface.  The  specific 
plasmonic  metasurface  shown  in  Fig.  la  consists  of  a  square  array 


of  asymmetric  unit  cells.  Each  unit  cell  consists  of  two  resonant 
gold  elements,  (i)  a  vertical  dipole  nanoantenna  and  (ii)  a 
horizontal  monopole  nanoantenna29’30,  connected  to  a  vertical 
plasmonic  wire  running  across  the  entire  metasurface.  In  the 
absence  of  the  near-held  interaction  between  the  nanoantennas, 
their  respective  plasmonic  resonances  appear  at  the  wavelengths 
2d=27 zc/cOy  and  2M=27ic/mx,  respectively.  The  near-held 
interaction  between  the  two  antennas  reduces  the  metasurface’s 
symmetry,  making  it  planar- chiral  because  of  the  absence  of 
mirror-rehection  symmetry14,3 1,32 .  The  minimal  mathematical 
description  of  the  nanoantennas’  interaction  with  each  other  and 
with  the  incident  light  wave  with  complex-valued  electric  held 
amplitudes  (Ex,  Ey)  is  given  by  the  following  set  of  equations 
obtained  from  the  temporal  coupled  mode  theory32-34: 


d  /flx\  Y  Ct>x  -K  ( Tx;rad2£x\ 

df  V  fly  /  \~K  ®y  J  \ay  /  \  zy-,rJdEy  ) 
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where,  ax>y  are  the  electric  dipole  moments’  amplitudes,  and 
«x,y=fflx.y-!(T-'+T-y;ohlJ  are  the  complex-valued 
frequencies  of  the  respective  plasmonic  resonances  that  take 
into  account  their  radiative  (Tx>y;rad)  and  Ohmic  (TXjy;0hm) 
lifetimes.  The  near-held  coupling  coefficient  k  between  the  two 
resonances  accounts  for  the  planar- chiral  nature  of  the 
metasurface,  giving  rise  to  the  proportional  CDT  (see  the 
‘Methods’  section  for  details): 
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Several  observations  can  be  made  from  equation  (2).  (a)  No  CDT 
is  possible  without  Ohmic  loss  in  at  least  one  of  the  two 
plasmonic  modes,  and  (b)  the  CDT  is  maximized  if  the  two 
modes  have  very  different  ratios  between  their  Ohmic  and 
radiative  lifetimes  (Tx;rad/Tx;0hm  ^  Ty;rad/Ty;0hm)  and  nearly  equal 
resonant  frequencies  (cox^coy).  The  latter  condition  informs  our 
specihc  choice  of  the  metasurface  because  the  monopole 
nanoantenna’s  mode  is  known21  to  have  a  much  longer 
radiative  lifetime  than  the  dipole  antenna’s  mode,  and  the  two 
can  be  spectrally  overlapped. 

To  conhrm  these  simple  analytic  predictions,  we  have  carried 
out  first-principles  COMSOL  simulations  of  the  metasurface 
shown  in  Fig.  1.  By  artificially  removing  Ohmic  losses  from  the 
simulation,  we  observed  from  Fig.  If  that  indeed  TR=TL  is  in 
accordance  with  equation  (1)  and  in  agreement  with  earlier 
theoretical  predictions13,19.  Therefore,  a  non-dissipative  planar 
metasurface  cannot  exhibit  optical  activity  despite  the  fact  that 
the  near-field  optical  intensities  produced  by  the  incident  LCP 
and  RCP  light  waves  are  markedly  different  in  their  spatial 
distribution  and  intensity,  as  can  be  observed  from  Fig.  lb-e. 
Ohmic  loss  drastically  changes  this  situation  and  induces 
considerable  CDT  as  shown  in  Fig.  lg.  The  qualitative  physical 
explanation  of  the  CDT  is  that  different  optical  intensity  and 
electric  current  distributions  induced  by  the  RCP  and  LCP  light, 
that  is,  weak  excitation  of  the  vertical  nanoantenna  for  the  LCP 
light  and  much  stronger  excitation  of  the  horizontal  nanoantenna 
by  the  RCP  light  as  shown  in  Fig.  lb,d,  leads  to  strong  CD  A 
(  AAcd  >  0)  in  the  spectral  proximity  of  the  plasmonic  resonance 
of  the  nanoantennas.  Therefore,  ATCD<0  is  expected  because  the 
transmission  of  the  RCP  light  is  weakened  by  the  light  absorption 
in  the  metal  (or,  as  it  would  be  relevant  for  microwave 
metasurfaces13,  in  the  substrate).  This  explanation  directly 
follows  from  the  analytic  model  given  by  equation  (1),  as 
detailed  in  the  ‘Methods’  section.  From  symmetry  considerations 
we  also  conclude  that  the  signs  of  A ACD  and  ATCD  are  reversed 
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Figure  1  |  Circular  dichroism  in  a  planar-chiral  plasmonic  metasurface  caused  by  Ohmic  losses,  (a)  Schematic  of  a  periodic  metasurface  (four  unit  cells 
shown)  illuminated  by  circularly  polarized  light  from  under  the  CaF2  substrate,  (b-e)  Surface  current  (arrows)  and  colour-coded  amplitude  (b,d)  and  grey¬ 
scale  phase  (c,e)  of  the  normal  electric  field  Ez  from  COMSOL  simulations  of  normally  incident  RCP  (b,c)  and  LCP  (d,e)  light  waves.  (f,g)  Transmission 
spectra  of  the  metasurface  for  the  cases  of  RCP  (blue  solid  line)  and  LCP  (red  dashed  line)  excitations  without  (f)  and  with  (g)  Ohmic  loss  in  the  metal. 
Simulation  parameters:  Px  =  Py  =  6  pm,  LM  =  1pm,  LD  =  3.3pm,  w=0.6pm,  g  =  0.4  pm.  Frequency-dependent  refractive  index  of  gold  n au  =  was 
obtained  from  ref.  44,  and  frequency-independent  refractive  index  ncaF2=L32  of  the  substrate  was  used. 


for  the  metasurface  enantiomer  obtained  by  the  y— >  —  y 
symmetry  transformation,  which  formally  reverses  the  sign  of  k 
in  equation  (1).  Note  that  the  effect  of  the  asymmetry  between  the 
substrate  (CaF2,  with  the  refractive  index  «caF2=l-32)  and  the  air 
superstrate  is  extremely  weak  as  observed  from  Fig.  If. 

Fabrication  and  optical  characterization.  Despite  its  apparent 
simplicity  and  qualitative  appeal,  this  conceptual  sequence — from 
different  near-field  distributions  produced  by  RCP/LCP  light  to 
differential  Ohmic  dissipation  to  CDT — has  never  been  experi¬ 
mentally  demonstrated  because  of  the  challenge  of  measuring  the 
near-field  distribution  and  Ohmic  loss  on  the  nanoscale.  Here,  we 
provide  the  experimental  evidence  on  the  dissipation’s  role 
in  inducing  chiral  dichroism  by  using  two  complementary 
AFM-based  techniques  for  directly  comparing  the  effects  of  the 
two  circular  light  polarizations  on  the  two  metasurface  enantio¬ 
mers:  (i)  transmission-mode  s-SNOM  for  measuring  the 
distribution  of  optical  fields  on  the  nanoantennas35-,  and 
(ii)  AFM-IR24-28  for  measuring  Ohmic  heating  of  the  metal. 

The  nanoscale  topographies  of  the  two  metasurface  enantio¬ 
mers  (referred  from  here  on  as  upper  and  lower)  are  shown  in 
Fig.  2a, b.  They  are  fabricated  on  top  of  a  mid-infrared  transparent 
CaF2  substrate  using  e-beam  lithography,  followed  by  gold 
deposition  and  the  lift-off.  Circularly  polarized  light  pulses  from  a 
tunable  quantum  cascade  laser  (QCL)  were  used  to  illuminate  the 
lower  enantiomer  (see  the  ‘Methods’  section  for  fabrication  and 
measurements’  details).  The  transmission  spectra  TR(2)  and  TL(2) 
shown  in  Fig.  2c  were  measured  across  the  8.25  jim<2<ll  pm 
spectral  range.  The  experimental  results  clearly  demonstrate 
strong  CDT,  which  is  most  pronounced  near  to  the  horizontal 
nanoantenna’s  resonance  XM  =  9  pm.  Specifically,  the  much  lower 
transmission  of  the  LCP  light  correlates  with  higher  Ohmic  losses 


of  the  LCP  for  the  lower  enantiomer.  Note  that  the  opposite  is 
true  for  the  upper  enantiomer,  see  Fig.  lg.  While  a  substrate  can, 
in  principle,  induce  small  CDT,  our  numerical  simulations  shown 
in  Fig.  If  indicate  that  the  effect  is  too  small  to  be  measurable. 

For  the  field- scanning  s-SNOM  experiments,  the  structure  was 
illuminated  through  the  CaF2  substrate  by  RCP  light  generated  by 
a  C02  laser  at  X  —  93  pm  (transmission-mode  s-SNOM35-37). 
Schematics  of  the  experimental  setup  are  shown  in  Fig.  3a  (see  the 
‘Methods’  section  for  measurement  details).  The  s-SNOM  is 
based  on  an  atomic  force  microscope  (AFM)  where  a  sharp 
dielectric  tip  locally  scatters  the  local  near  fields  on  the  sample 
surface.  While  the  sample  is  scanned,  the  scattered  light  is 
detected  interferometrically  using  a  p -polarized  reference  beam, 
and  amplitude  and  phase  maps,  Ep  and  <pp,  of  the  near-field 
distribution  of  both  enantiomers  are  obtained.  In  case  of  the 
presented  structure,  these  maps  mainly  correspond  to  the  vertical 
near-field  component,  Ez  (see  Fig.  5  in  the  ‘Methods’  section 
section).  The  experimental  results  presented  in  Fig.  3b, c  clearly 
reveal  that  the  near-field  optical  intensity  distribution  strongly 
depends  on  the  handedness  of  the  enantiomer  for  a  given 
handedness  of  incident  light.  In  more  detail,  we  observe  the 
fundamental  dipolar  mode  on  the  dipole  antenna  for  both 
enantiomers,  exhibiting  strong  amplitude  signals  at  the  rod  ends 
and  a  phase  jump  at  the  rod  center.  The  monopole  antenna  shows 
strong  amplitude  signal  only  for  the  upper  enantiomer,  while  the 
fields  completely  vanish  at  the  monopole  antenna  for  the  lower 
enantiomer.  From  the  abovementioned  symmetry  considerations, 
it  follows  that  the  response  of  the  lower  enantiomer  to  the  RCP 
light  is  identical  to  the  response  of  the  upper  enantiomer  to  the 
LCP  light.  We,  thus,  have  experimentally  demonstrated  for 
the  first  time  that  the  near-field  intensity  distribution  near  a 
given  planar- chiral  metasurface  strongly  depends  on  light’s 
handedness32,38. 
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Figure  2  |  Far-field  measurements  of  CDT  in  a  planar-chiral  metasurface.  (a,b)  AFM  image  of  the  enantiomers  fabricated  on  top  of  CaF2  substrate.  The 
two  enantiomers  (top  and  bottom)  are  mirror  images  of  each  other,  (c)  Experimental  far-field  transmission  for  the  structure  shown  in  the  inset  (and  also 
corresponding  to  (b))  revealing  the  CDT  caused  by  the  Ohmic  loss.  Open  circles  represent  the  experimental  points  for  LCP  (blue)  and  RCP  (red) 
excitations;  solid  lines  are  for  guiding  the  eye.  Inset  shows  SEM  image  of  the  enantiomer  studied. 


Figure  3  |  s-SNOM  measurements  of  the  optical  field  above  two  planar-chiral  metasurfaces  illuminated  by  RCP  light,  (a)  Schematics  of  the 
experimental  setup,  (b)  Measured  amplitude  |EP  |  and  (c)  phase  of  the  normal  electric  field  Ep  at  X  =  9.3  pm.  White  dashed  lines  in  b,c  outline  the  geometry 
of  the  structure.  BS,  beam  splitter;  Det,  pseudoeheterodyne  detection;  P,  linear  polarizer;  PM,  parabolic  mirror;  PRM,  infrared  phase  retarding  mirror. 


The  qualitative  explanation  of  this  effect  rests  on  the 
observation  that,  depending  on  the  enantiomer,  the  vertical 
dipole  nanoantenna  (which  is  primarily  excited  by  the  Ey 
component  of  the  incident  laser  beam)  either  polarizes  or 
de-polarizes  the  horizontal  monopole  nanoantenna  excited  by 
the  £x  field  component.  This  effect  is  most  pronounced  for  CP 
light  when  the  frequencies  of  the  two  modes  are  spectrally 
matched  to  each  other  and  to  the  incident  laser’s  frequency 
(co^cox^coy),  so  that  the  phases  of  ax>y  with  respect  to  £Xjy  are 
approximately  the  same.  Therefore,  the  ±  7i/2  phase  shift  between 
£x  and  Ey  (depending  on  light’s  handedness)  results  in  either 
constructive  or  destructive  interference  depending  on  the 
enantiomer  (that  is,  the  sign  of  k  in  equation  (1)).  This 
theoretical  argument  is  validated  by  the  excellent  agreement 
between  the  experimental  results  presented  in  Fig.  3b, c,  and  the 
results  of  numerical  modelling  presented  in  Fig.  lb-e. 


Measurements  of  polarization-dependent  Ohmic  heating.  We 

proceeded  to  demonstrate  that  the  distinct  optical  intensity 
distributions  produced  by  the  two  circular  light  polarizations 
shown  in  Fig.  3  give  rise  to  correspondingly  different  Ohmic 
losses  inside  the  plasmonic  metasurface.  This  was  accomplished 
by  using  the  AFM-IR  technique,  which  relies  on  detecting 
thermal  expansion  of  polymers  and  other  materials  on  nanos¬ 
cale  by  observing  cantilever  deflection  of  the  AFM.  While  AFM- 
IR  technique  has  originally  been  developed  for  nanospectro¬ 
scopy  of  organic  compounds24,25,28,  it  can  also  be  used  to 
image  local  heat  generation  in  nanoantennas  by  observing  the 
expansion  of  the  gold  itself,  or  an  infrared-transparent  polymer 
film  in  contact  with  nanoantennas26,27.  Polyethylene  has 
virtually  no  infrared  absorption  at  the  laser  wavelength,  and 
its  heating  is  produced  entirely  by  Ohmic  dissipation  in 
nanoantennas. 
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Figure  4  |  AFM-IR  measurements  of  circularly  dichroic  thermal  expansion  of  the  metasurface,  (a)  Experimental  setup.  Light  from  a  QCL  passes  through 
a  quarter-wave  plate  (QWP)  and  illuminates  the  sample  from  CaF2  substrate.  Local  Ohmic  heating  of  nanoantennas  causes  local  heating  and  expansion  of 
the  polyethylene  (PE)  film  and  leads  to  cantilever  deflection,  which  is  detected  by  observing  laser  tracer  beam  deflection  from  cantilever,  using  position- 
sensitive  photodetector  (PSPD)  and  a  lock-in  amplifier,  referenced  by  the  QCL  pulse  repetition  frequency,  (b)  COMSOL  simulations  of  the  temperature 
increase  distribution  in  the  polyethylene  film  at  the  end  of  a  square  QCL  pulse  with  time  duration  T  =  200 ns  and  peak  intensity  /  =  1  kWcm~2  tuned  to 
2  =  9.1  pm,  corresponding  to  the  experimental  conditions,  (c)  AFM-IR  cantilever  deflection  on  top  of  a  PE-coated  sample  excited  with  RCP  (left  panel)  and 
LCP  (right  panel)  laser  pulses  at  normal  incidence  through  the  CaF2  substrate. 


For  our  measurements,  we  evaporated  a  100-nm  thick  layer  of 
polyethylene  on  top  of  the  metasurface  made  of  the  lower 
enantiomer  shown  in  Fig.  2b.  The  layer  of  polyethylene  helps  us 
to  improve  the  photo-expansion  signal  because  the  linear  thermal 
expansion  coefficients  of  polyethylene  are  over  an  order  of 
magnitude  higher  than  that  of  gold.  Circularly  polarized  200-ns- 
long-light  pulses  at  2  =  9.1  pm  from  a  tunable  QCL  were  used  to 
illuminate  the  sample  through  the  CaF2  substrate  at  normal 
incidence  with  1=1  kW  cm -  2  peak  intensity  on  the  sample 
(see  Fig.  4a  and  the  ‘Methods’  section  for  the  details  of  the 
experimental  setup  and  measurements).  Figure  4b  shows  the 
results  of  COMSOL  simulations  of  the  temperature  distribution 
in  the  polyethylene  film  at  the  end  of  the  LCP  and  RCP  laser 
pulses.  Because  thermal  diffusion  length  in  polyethylene  is  well 
below  100  nm  during  the  laser  pulse  polyethylene,  heating  is 
found  to  be  highly  localized  and  to  closely  follow  the  distribution 
of  Ohmic  dissipation  (that  is,  the  current  density)  in  the  metal 
nanoantennas  shown  in  Fig.  lb-e. 

Figure  4c  shows  the  experimentally  measured  AFM  cantilever 
deflection  amplitude  at  different  areas  of  the  sample.  The 
cantilever  deflection  is  directly  proportional  to  temperature 
increase  in  the  sample  during  the  laser  pulse.  The  experimental 
data  are  in  excellent  agreement  with  theoretical  predictions 
shown  in  Fig.  4b.  It  confirms  that  the  magnitude  and  spatial 
distribution  of  the  Ohmic  heating  of  a  chiral  2D  metasurface 
markedly  depends  on  the  handedness  of  light.  For  example,  the 
vertical  dipole  nanoantenna  is  extinguished  for  resonant  LCP 
illumination.  Considerably  less  polarization-induced  asymmetry 
is  induced  by  off-resonant  illumination  (see  Supplementary  Fig.  1 
for  thermal  expansion  images,  and  Supplementary  Note  1  with 
Supplementary  Fig.  2  for  theoretical  interpretation).  These  results 
represent  the  first  direct  experimental  evidence  of  Ohmic  loss 
origin  of  CDT  in  a  planar- chiral  plasmonic  metasurface. 

Discussion 

In  conclusion,  by  using  three  distinct  experimental  techniques 
(two  near-field  and  one  far-held),  we  have  demonstrated  that 
optical  activity  in  planar- chiral  plasmonic  metasurfaces  is 


unambiguously  correlated  with  circularly  dichroic  Ohmic  loss, 
which  is  in  turn  caused  by  circularly  dichroic  near- held 
distribution.  The  direct  mapping  of  nanoscale  optical  helds  and 
absorption  rates  were  carried  out  for  two  chiral  metasurface 
enantiomers  illuminated  by  circularly  polarized  infrared  light. 
Although  the  importance  of  this  work  is  primarily  fundamental,  it 
also  paves  the  way  for  practical  applications,  such  as  the 
development  of  novel  optical  detectors  for  Stokes  parameters’ 
polarimetry,  as  well  as  for  surface-enhanced  enantiomeric  sensing 
of  chiral  biological  molecules3 1,38>39.  The  design  principles  behind 
achieving  the  large  circularly  dichroic  effects  in  plasmonic 
metasurfaces  that  are  established  and  experimentally  verihed  in 
this  study  will  be  of  great  interest  to  a  variety  of  researchers, 
especially  in  the  held  of  biological  sensing  and  novel  infrared 
optical  devices. 

Methods 

Fabrication  of  plasmonic  planar-chiral  metasurfaces.  Planar-chiral  plasmonic 
metasurfaces  were  fabricated  on  infrared-transparent  CaF2  substrates  using 
electron  beam  lithography,  followed  by  a  metal  lift-off.  Polymethyl  methacrylate 
electron  beam  resist  spun  at  1,700  r.p.m.  for  30  s  was  pre-baked  for  1  min,  followed 
by  the  deposition  of  a  thin  ( ~  3-5  nm)  gold-paladium  conductive  layer  using 
sputtering.  The  desired  structures  were  written  using  a  JEOL  6000  EBL  system 
at  450nCcm-2,  and  developed  in  1:3  MIBK:IPA  solution  for  60s.  After  devel¬ 
opment,  the  sample  was  washed  with  IPA  and  water,  followed  by  a  brief  oxygen- 
plasma  cleaning  step  (100  W  for  5  s)  to  prepare  for  metal  deposition.  A  60-nm 
thick  gold  layer  was  thermally  deposited  following  a  5-nm  thick  Cr  adhesion  layer. 
The  sample  was  then  immersed  in  acetone  for  lift-off,  followed  by  oxygen-plasma 
cleaning. 

s-SNOM  measurements  of  the  field  distribution  at  the  metasurface.  Our 

s-SNOM  setup  (see  Fig.  3a  for  a  schematic)  is  based  on  an  atomic  force  microscope, 
where  commercial  silicon  tips  (NanoWord,  Arrow-NCR-50)  were  used  to  locally 
scatter  the  near  fields  on  the  sample  surface.  The  sample  and  tip  were  illuminated 
from  below  at  normal  incidence  with  a  weakly  focused,  right-handed  circularly 
polarized  C02  laser  beam  at  a  wavelength  2  =  9.3  pm  (transmission-mode 
s-SNOM35-37),  thus  ensuring  homogeneous  illumination  of  a  relatively  large  area 
on  the  sample  without  phase-retardation  effects.  Using  a  parabolic  mirror,  we 
collected  the  scattered  light  in  the  x-direction  at  an  angle  of  60°  from  the  surface 
normal.  The  scattered  light  was  then  interfered  with  a  phase-modulated,  vertically 
(z-)  polarized  reference  beam  at  the  infrared  detector  and  recorded  simultaneously 
with  the  sample  topography.  Background  contributions  could  be  fully  suppressed 
by  vertical  tip  oscillation  at  a  frequency  Q  =  250  kHz  (tapping-mode  AFM)  and 
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Figure  5  |  Comparison  between  numerical  calculations.  Displaying  Exr  Ey,  and  Ez  near-field  components  at  a  height  of  50  nm  above  the  structures  (three 
left  subplots)  and  the  experimentally  recorded  near-field  maps,  showing  the  p-component  of  the  scattered  field,  Ep  (right  subplots).  The  comparison  shows 
that  the  experimental  maps  mainly  yield  the  vertical  near-field  component,  Ez.  (a,b)  correspond  to  the  two  enantiomers  of  the  metasurface.  The  colour 
scheme  for  the  simulations  results  is  the  same  as  for  the  experimental  results. 


subsequent  higher  harmonic  demodulation  of  the  detector  signal  at  3(2  (ref.  40).  In 
this  case,  the  vertical  polarization  of  the  reference  beam  selected  the  p-component 
of  the  scattered  light.  Using  a  pseudoheterodyne  detection  module 
(www.neaspec.com)  the  amplitude  |£p|  and  phase  cpp  of  the  scattered  light  were 
measured  for  each  scanning  point. 

Recently,  it  was  found  that  there  are  different  scattering  mechanism  in  s-SNOM 
imaging  of  metal  antennas:  (i)  direct  scattering  of  the  antenna  near  fields  by  the  tip 
into  the  far  field  (detector)  and  (ii)  scattering  of  the  antenna  near  fields  by  the  tip 
via  the  antenna  itself41,42.  When  imaging  resonant  linear  antennas,  excitation  with 
s-polarization  and  detecting  the  s-component  of  the  scattered  light  is  mainly  based 
on  the  mechanism  (ii),  while  detecting  the  p-component  is  mainly  based  on  the 
mechanism  (i),  where  the  latter  was  found  to  essentially  yield  the  z-component,  Ez, 
of  the  near  fields  on  the  antenna  surface42.  In  case  of  the  presented  metasurface,  the 
p-component  of  the  scattered  field  mainly  yields  the  vertical  (z)  component,  Ez,  of 
the  near  fields,  as  it  was  found  by  comparing  the  experimental  images  with 
numerical  calculations37.  The  vertical  (z)  component,  Ez,  is  closely  related  to  the 
magnitude  and  phase  of  the  surface-charge  density  in  metal  antennas,  so  that  the 
near-field  maps  shown  in  Fig.  3  allow  for  unambiguous  identification  of  the 
antenna  modes. 


Comparison  of  numerical  field  profiles  and  s-SNOM  measurements.  As  the 

s-SNOM  imaging  relies  on  interferometric  detection  of  the  scattered  light  using  a 
p-polarized  reference  beam,  the  correspondence  between  the  experimental  mea¬ 
sured  maps  of  the  field  amplitude,  Ep  and  phase,  (pp,  and  the  numerically  calculated 


fully  vectorial  field  profiles  should  be  established.  As  can  be  seen  from  Fig.  5,  in 
case  of  the  structure  presented,  comparison  of  the  experimental  and  numerical 
maps  clearly  reveals  that  the  measured  field  components  correspond  to  the  vertical 
Ez  near-field  component  of  the  electric  field. 


AFM  measurements  of  thermal  expansion.  A  pulsed  mid-infrared  QCL  (Day¬ 
light  Solutions,  Inc.;  tuning  range  of  900-1,200  cm  -  x)  operating  at  9.1-|im 
wavelength  was  used  to  illuminate  the  sample.  The  light  from  the  QCL,  originally 
linearly  polarized,  was  converted  to  LCP  or  RCP  light  using  an  achromatic  quarter- 
wave  plate  (Altechna  Inc.,  2-IRPW-ZO-L/4-8000-C).  The  laser  was  operated  with 
200-ns  pulses  with  the  peak  power  of  300  mW  at  a  repetition  frequency  of 
~  180  kHz,  in  resonances  with  the  second  bending  mode  of  the  AFM  cantilever.  To 
increase  the  photo-expansion  signal,  a  100-nm  thick  polyethylene  film  was  ther¬ 
mally  evaporated  on  top  of  the  nanoantennas.  The  sample  was  positioned  on  the 
AFM  stage  with  the  polyethylene-coated  metal  nanoantennas  being  on  top,  right 
below  the  AFM  tip.  Laser  radiation  was  incident  normally  onto  the  sample  from 
below  from  the  CaF2  substrate.  The  beam  was  focused  to  a  100-jim  radius  spot 
using  ZnSe  lens  with  the  convergence  half-angle  of  4°.  The  QCL  peak  intensity  on 
the  sample  surface  was  estimated  to  be  ~  1  kW cm-2. 

The  AFM  system  was  operated  in  contact  mode.  Gold-coated  ContGB-G 
(Budget  Sensor)  AFM  cantilevers  were  used  for  measurements.  Thermal  sample 
expansion  during  the  laser  pulse  exerted  a  force  on  the  AFM  cantilever  tip  that 
leads  to  cantilever  deflection.  The  cantilever  deflection  amplitude  was  measured  by 
a  position-sensitive  photodetector  and  sent  to  a  lock-in  amplifier  (Stanford 


6 


NATURE  COMMUNICATIONS  |  7:12045  |  DOI:  10.1038/n com msl 2045  |  www.nature.com/naturecommunications 


NATURE  COMMUNICATIONS  |  DPI:  10.1038/ncomms12045 


ARTICLE 


Research,  SR844)  using  QCL  pulse  repetition  rate  as  a  reference.  The  lock-in 
integration  time  was  set  at  10ms.  The  output  from  the  lock-in  amplifier  was  used 
to  form  64  x  64  pixel  images  shown  in  Fig.  3b  by  raster  scanning  the  sample  at  the 
speed  of  0.2  Hz.  In  our  Multiphysics  COMSOL  simulations,  which  combined  the 
simulation  of  the  laser-induced  metasurface  heating  and  thermal  conduction,  the 
following  thermal  diffusivities  were  used  for  gold,  CaF2  substrate  and  polyethylene: 
DAu  =  1.28  x  10_4m2s_  1,  Dsub  =  3.6  x  10_6m2s_  and  DPE  =  10_7m2s_  \ 
respectively. 

We  have  also  measured  (not  shown)  the  metal  expansion  directly,  that  is, 
without  the  PE  polymer  layer,  at  the  resonant  (/R  =  9.1  pm)  and  at  the  non¬ 
resonant  (2nr=  10.3  pm)  wavelengths.  Strong  handedness-dependent  expansion  of 
the  monopole  antenna  and  the  adjacent  area  of  the  continuous  nanowire  was  found 
only  at  X  =  /R,  whereas  much  weaker  metal  expansion  was  observed  at  X  =  zNR  for 
both  circular  light  polarizations. 


corresponds  to  the  matched  plasmonic  resonances  of  the  two  nanoantennas 
(cox  =  coy=coR)  and  the  following  hierarchy  of  radiative  and  non-radiative  losses: 
Ay; rad  ^  T*,ohm;  ^y;Ohm  <  Tx;rad.  In  this  limiting  case  the  expression  for  the 
absorption  rates  at  resonance  ( co  =  coR)  can  be  simplified  by  noting  that  most  of  the 
absorption  happens  at  the  monopole  antenna: 


Tx;Ohm 

1  K--\/^y;rad^x;rad 

^x;rad 

1  +  JC2Ty;radTx;0hm 

Tx;Ohm 

1  +  ^-y/^y;rad^x;rad 

G:;rad 

1  +  K:2Ty;radTx;0hm 

It  clearly  follows  from  these  expressions  that  (a)  larger  loss  (and,  corres¬ 
pondingly,  lower  transmission)  occurs  for  the  RCP  light  in  the  case  of  the  upper 
enantiomer  with  k>0  shown  in  Fig.  la;  (b)  the  situation  reverses  for  the  lower 
enantiomer  (see  Fig.  2,  bottom  row)  which  is  characterized  by  k  <  0  as  explained  in 
the  text. 


Far-field  measurements  of  circular  dichroism  in  transmission.  A  pulsed  mid- 
infrared  QCL  (Daylight  Solutions,  Inc.;  tuning  range  of  900- 1,200  cm  -  x)  was  used 
as  a  light  source.  Linearly  polarized  light  from  QCL  passed  through  achromatic 
quarter-wave  plate  (Altechna  Inc.,  2-IRPW-ZO-L/4-8000-C)  for  polarization 
control  and  a  numerical  aperture  0.5  reflective  objective  lens  (Edmund  Optics, 
Adjustable  ReflX  Objective  36X/0.5NA  IR  FIN)  to  the  sample.  Transmitted  light 
from  the  sample  was  collected  and  collimated  by  a  3-in  focal  length  plano-convex 
ZnSe  lens  and  refocused  by  another  identical  lens  on  a  liquid  nitrogen-cooled 
mercury-cadmium-telluride  photodetector.  The  transmission  spectral  data  from  the 
metasurface  were  normalized  using  transmitted  signal  from  a  bare  CaF2  substrate. 

Analytic  coupled  mode  theory  of  circular  dichroism  in  planar-chiral  meta¬ 
surfaces.  The  temporal  coupled  mode  theory  expressed  by  equation  (1)  is  the 
minimal  description  of  a  metasurface  in  terms  of  its  electromagnetic  resonances 
that  enables  the  calculation  of  transmission  and  reflection  coefficients  for  arbi¬ 
trarily  polarized  normally  incident  light  wave  with  electric  field  amplitudes  (£x,  Ey). 
For  example,  the  reflectivity  matrix  R  is  obtained  from  the  following  equation: 


by  continuity,  the  transmission  matrix  is  given  by  T  =  I  +  R  ,  where,  I  is  a  unity 
matrix.  The  state  of  polarization  of  the  LCP/RCP  light  propagating  in  the 
z-direction  is  given  by  the  following  vectors: 


The  CDT  given  by  A  TCD=  TR—  TL  can  be  expressed  as 
A TCD =2Im { T*x Txy  —  T^T^},  revealing  that  if  the  transmission  matrix  can  be 
diagonalized  by  rotating  its  axis  by  an  arbitrary  angle  9,  then  ATCD  =  0.  The 
frequency-dependent  rotation  angle  9  is  given  by  the  following  expression: 


1  _  G;rad(&A  (&)  T  ^y;rad(^  &^y)  "f"  *'(^x;rad/^x;Ohm  ^y;rad/^y;Ohm) 

tan  29  2K^/TX;rad^y;rad 


(5) 


from  which,  it  follows  that  a  real- valued  9  exists  as  long  as  T^rL^yi^’  that  is, 
the  Ohmic  losses  vanish14,19. 

The  vanishing  of  the  CDT  in  the  absence  of  Ohmic  losses  can  be  understood  by 
expressing  ATCD  =  TLR—  TRL,  where,  TLR  and  TRL  are  the  transmission  conversion 
efficiencies  from  RCP  to  LCP,  and  from  LCP  to  RCP,  respectively13,43.  The 
markedly  different  charge  distributions  on  the  two  nanoantennas  under  RCP  and 
LCP  illuminations  (see  Fig.  lb-e  for  amplitude  and  phase  information)  are 
formally  equivalent  to  |ajcp|  produced  under  resonant  RCP  illumination  being 
much  larger  than  |ayCP|  produced  under  either  LCP  illumination  (note,  from 
Fig.  lb  that  | ayCP  |  is  negligible  under  RCP  illumination)  .  Nevertheless,  the 
x-polarized  monopole  antenna  resonance  with  the  amplitude  |flxCP|  couples  to 
LCP  light  with  the  same  efficiency  as  the  y-polarized  dipole  antenna  resonance  with 
the  amplitude  |  4CPi «  |4T  couples  to  RCP  light.  The  reason  for  which  is  the 
much  lower  radiative  rate  of  the  monopole  antenna:  1  /  Tx;rad  -Cl/ Ty;rad-  This 
qualitative  argument  explains  why  TLR=  TRL  and,  therefore,  ATCD  =  0  despite 
strong  differences  in  the  near-field  distributions  is  shown  in  Fig.  1. 

For  non-vanishing  Ohmic  losses,  the  total  absorption  coefficients  for  the  LCP 
and  RCP  illuminations  are 


Al 


2 


Tx;Ohm  \k2  —  (co  ~  COx)  ( CO  —  COy)  |2 


\2 


Ty;Ohm  |  K2  —  (CO  —  COx)  (CO  —  COy)  | 2 


(6) 


ar  = 


Tx;Ohm  \k2  —  {co  —  COx)  (co  —  C0y)  |  Ty;Ohm  \k2  —  (a)  —  COx)  (co  —  COy)  | 


(7) 


respectively,  where  the  two  terms  in  the  sum  represent  the  contributions  of  the 
monopole  and  dipole  nanoantennas,  respectively.  The  simplest  limiting  case 


Data  availability.  The  data  that  support  the  findings  of  this  study  are  available 
from  the  corresponding  author  on  request. 
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